On surface energy of graphene and carbon nanoclusters 
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Theoretical study of graphite (graphene) edge is done. The most stable edge orientation is 
calculated to be a zigzag [110] edge. Possible applications of the result to the formation of different 
graphitic structures are discussed. 



Q ' I. INTRODUCTION 

Recently, attention of physicists and chemists was recalled to the problem of reconstruction of edges of planar 
(and curvilinear) graphite. Considerable amount of articles on nanocarbons |l[, surmounted with recent results on 
nanotubes || and graphite polyhedral crystals j|, shows that this question is far from being clear, despite intensive 
experimental and theoretical research has been conducted. 

The aim of this paper is to get theoretical understanding of the formation of facets and edges of graphite, graphene 
(monolayer of graphite) and small carbon sheet fragments. I will use simple analytical approach implementing the 
i-rt ! continuum energetics model H assumptions. This approach bases on the count of dangling bonds along the free edge 
£h ■ of graphene. That gives an energy of the facet[] as a function of its orientation. The hexagonal lattice of the graphene 
has a three-fold symmetry and allows me to consider only a sector of the full plane between and 2-7r/6. 

The function of specific facet energy, F(uo) = E(u))/L, is a periodic function of the orientation angle, ui. At the 
fixed u) it has a weak dependence on the facet length, L, (in 3D case it could depend on the length of the sheet edge 
or on the free surface area, which is defined by the specific type of misorientation). This additional dependence is 
fully neglected if considering the continuous edge, owing to full neglect of the atomic structure. A microscopic view 
is that at the scale of interatomic distance, b, one can not choose an arbitrary orientation of the edge. The natural 
(f) | approach to the description of free misoriented facet is to consider a number of one-atom kinks distributed along the 
C — ■ edge so that their density depends on the (continuous) orientation angle. 

The orientational dependence of the facet energy appears to be important for the edge zipping model of the single- 
wall nanotube nucleation M. In order to understand the experimental fact that synthesized SWNTs had indices 
[n,n±(m] close to armchair helicity n the zipping of the graphite edges with different orientations has to be considered. 



II. CHIRALITY DEPENDENT FORMATION ENERGY 



The surface energy of free graphene facet depends on its orientation via the energy of dangling bonds (DBs). In 
this context, the term free facet means that there is no special reconstruction of the edge bonds. Of course, the 
reconstruction can change the facet energy substantially. I consider some type of the 3D facet reconstruction in Ref. 

. £h O • The reconstruction of the 2D facet is somewhat more complicated as the system possesses less degree of freedom. 

S^ • It will be discussed at length elsewhere. 

The microscopic picture of the free edge is presented in Fig.Rl (upper). Along the graphene edge the dangling bonds 
are schematically drawn. The orientation of the main edge is chosen along [210] zigzag direction PI which is defined 
respectively to the basis lattice unit vectors c\ and oz (the third unit vector S3 is normal to the graphite plane, so the 
corresponding lattice number is always equal to for graphene) . The density of DBs increases at such discontinuities 
of the edge as the kink shown in the Figurcfj]. 



* E-mail: rotkin@theory.ioffe.rssi.ru 

1 Henceforth I will use the term facet for the 2D planar face of 3D graphite crystal and ID straight edge of 2D graphene sheet. 
That will not lead to misunderstanding because the latter is a particular case of the former when the number of layers equals 



one 

2 



To mark a SWNT I follow the conventional index agreement proposed in Ref. M. 
3 There are multiple possible choice of the unit vector diads; the second (tilde) set, shown in the Figure 1 on the right, gives 
the edge symbol [110]. Both choices are equivalent to [110] 
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FIG. 1. Misoriented edge of the graphene layer. The exact direction of the 2D-crystal-face, A, is defined as an envelope 
along the kinked edge. Any edge, being different from zigzag and armchair basal types, is, at the atomic scale, a series of kinks 
and steps of these basal faces. Upper: Details of the single kink at the zigzag edge and the unit cell geometry of the graphene. 



The total energy of DBs constitutes the free surface energy: 

-Sedge = -/Viayers v{uj) L Eb, 



(1) 



here the edge density of DBs, i'(w), the edge length, L, and the number of layers, -/Vi ayors , give the total number of 
DBs. Eb ~ 2.355 eV is the single DB energy fl. I do not consider here any passivation/contamination of the graphite 
surface, which could change Eb- 

The less the edge density, the more stable the surface. It is the zigzag face normal to [110] (or symmetry equivalent) 
direction fllPf . The edge, being scrolled and zipped (as shown in Fig.gJ), will form a cylindrical sleeve (SWNT nucleus). 
Its axis is along direction [110] in the graphene plane. This corresponds exactly to armchair nanotube conformation. 
I speculate that it gives a plausible explanation of the preferable formation of armchair SWNTs at certain conditions. 
The direct observation of graphene scrolling was published in Ref. |llHl3j| . It gives also an experimental support that 
the stable graphene edge is the zigzag one. 




FIG. 2. Simulation of the graphene edge zipping (after [14] 



The edges of zigzag and armchair types have the minimum and maximum densities of the DBs. It is possible to 
expand the DB density around extremum points and write an interpolation formula in between. Surprising cuspn is 
found for the minimum, where the density changes linearly in to, the chirality angle. The latter is a standard measure 
of how far does the SWNT axis deviate from the high symmetry direction, for example, the zigzag direction. Because 
my edge structures have ID geometry close to tubular one (Fig. @), it is natural to characterize their misorientation 
with the same chirality angle. 

Let me define the edge density through the nanotubc-likc indices n,m which are the numbers of unit vectors 
constituting the vector A joining the left and right ends of the considered edge structure [n m 0] (Fig.0). The DB 
density is the number of boundary sites divided by the length of the boundary: 



\A\ = Vsby/rf' 



ms 



(2) 



4 A similar nonmonothonic behavior is known for the function of the surface energy of perfect crystals for a century at least 
(G.S.Wulff, Z.Krist.Mineral. v. 34, 449, 1901). 



where b — \c\/y/3 ~ 1.41 A is the carbon bond length. The infinitesimal small shift of the vector A =[n,n] (or [n,0]) 
defining the armchair (or zigzag) basal edge results in the change of the edge length and the number of DBs. I expand 
both in series on the small parameter 1/n ex oj. For an armchair edge it reads as: 

2n+l _2_ / J_ l_ J_ _3_ J_ \_ 2_ ( J_^ 2 \ 

VA ~ 3&n^l + l/n + l/3n 2 ~ 36 V. 2n + 2n 6n 2 + 8n 2 4n 2 + ' ' 7 36 V ^ '")' U 



Naturally, the first order terms cancel while the second order terms give the square dependence in the chirality angle, 
oj ex — , around the maximum. The coefficient is fc^ = — - 
The non-zero expansion terms appear in the first order: 



oj ex i, around the maximum. The coefficient is k^ — — i- The behavior of a nearly zigzag edge density is different. 



71+ 1 1/11 \ 1 / 1 \ 

v z = —= =- H ... = —=- 1 H =w + ...). (4) 

VSfenv/i + V^+l/" 2 V36 V n In ) 736 V 273 / 

Hence, the linear dependence coefficient is sz = H — y= and the density of DBs has a cusp near the point of the 
minimum. 

An edge, which has the intermediate crystallography index [n m 0] is nothing but a series of kinks (single atom 
steps) along the basal edge of a definite type (FigjI]). 

Let me consider the (approximately) circular fragment of graphene, which is naively expected to deliver the minimum 
to the surface energy in a continuum approximation. The actual surface energy will oscillate along the perimeter with 
the period of ir/3 between maxima corresponding to equivalent armchair edges (Figure |3|). 
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FIG. 3. Dangling bond density of the graphene edge vs. chirality angle. The surface energy of the graphene fragment is 
periodic function of u> because of v[u). Dotted lines represent the approximations at the extremum points given by Eqs.(3,4). 

At zero temperature, the purely zigzag facing of a 2D-crystal is expected. At the thermodynamic equilibrium, the 
zigzag-like edges are formed and the energetic cost for the vicinal edge is proportional to exp(— i/LEb/T), where v 
depends on the orientation (facing angle) linearly, see Eq.(Q). 

Let me quantify this. The energy gain for the nucleation of the vicinal zigzag edge [n, <5n, 0] is as follows: 

8E(u>) = Mayors 5v(w) E b L ~ Mayers L V Z E b S Z LO ~ ^^ E b . (5) 

Here I subtract the DB energy of the optimum zigzag edge, Mayers Lvz E b . The chirality equals oj = Sn/n. 

For the zipping model of SWNT nucleation one has Mayers = 2, it yields the following exponentially small probability 
of the nucleation of chiral SWNT with indices [n, n±<5n] at the temperature T: 
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w D exp[- — — £*] = w e UJ c i (6) 

where w is the probability of nucleation of armchair [n,n] SWNT, oj c = — — is a maximum chirality at given L 

L E b 
and T. The higher the temperature, the shorter the nucleus length. Thus the oj c has a maximum as a function of 

temperature. Let me estimate L ~ 50A at T ~ 2000 K. Then the oj c ~ 0.006 and the probability to nucleate [10,11] 

tube is 10" 7 smaller than for [10,10] SWNT 

Let me consider the edge closest to zigzag one. The smallest edge defect is a couple of kinks (per total nucleus 

length) with the energy equals E b . This minimum energy gain corresponds to the relative nucleation probability 

w/w ~ e- Eb ' T ~ 10~ 6 at T = 2000K. 



The calculated probability of chiral SWNT nucleation is too small to describe chirality scattering of the experiment. 
The real growing edge shape is likely far from equilibrium (ideal straight line). Therefore, the kinetics of graphene 
edge growth defines both the length and the chirality of the nucleus. I did not consider here the dangling bond 
passivation. Whence the kink energy is Eh/2, the passivation increases the maximum allowed chirality that could also 
explain the experimental scattering of SWNT angles. 

III. SUMMARY 

The microscopic model describing the energy of the edge formation for different types of graphitic structures is 
developed. The most stable edge has a zigzag orientation owing to the smallest density of the dangling bonds along 
the facet. The DB density is calculated depending on the misorientation angle for the facet of the arbitrary direction. 

In respect to proposed earlier SWNT nucleation mechanism via formation of a cylindrical sleeve along the graphene 
edge, this paper allows to speculate that the prevalent tubes have almost the armchair type which is complementary 
to the most stable zigzag graphene edge. The prevalent radius of such a sleeve nucleus is close to the optimum 
radius R m of nano-arch structure considered in 0] . My estimation gives a lower limit of the uniformity of the SWNT 
synthesis yield. In ideal conditions, the novel nucleation mechanism predicts formation of [11,11] SWNT (and no 
chirality scattering is expected) which is best approximation of continual energetics model to the most abundant 
[10,10] SWNT observed experimentally. 

Note, that in this paper I did not consider any kinetics of the formation, that might change the result. 
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